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Abstract 

The properties of a particle diffusing on a one-dimensional lattice where at 
each site a random barrier and a random trap act simultaneously on the par- 
ticle are investigated by numerical and analytical techniques. The combined 
effect of disorder and traps yields a decreasing survival probability with broad 
distribution (log- normal) . Exact enumerations, effective-medium approxima- 
tion and spectral analysis are employed. This one-dimensional model shows 
rather rich behaviours which were previously believed to exist only in higher 
dimensionality. The possibility of a trapping-dominated super universal class 
is suggested. 
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I. INTRODUCTION 



One-dimensional models are widely used in the physics of disordered systems This 
is because on one hand they are often good representatives of higher dimensional models 
and on the other hand they are much easier to handle. In the last decades the importance 
of investigating simple diffusion in the presence of disorder and trapping has been widely 
appreciated as toy model for complex systems including migrations of optical excitations 
0, polymer physics and diffusion-limited binary reactions J|. See e.g. [H and || for an 
exhaustive review. 

From the mathematical viewpoint the analytical solution of the diffusion problem of 
particles diffusing before getting completely absorbed by permanent traps, has been proven 
an extremely difficult task and only an asymptotic solution in presence of uncorrelated 
disorder for the survival probability (to be defined below) could be given using sophisticated 
techniques 0. 

A variation of this problem in the presence of strongly correlated (percolating) disorder, 
was also recently numerically investigated |||| and it was observed that whenever the total 
survival probability is an erratically decreasing function (i.e. detailed balance is violated) 
new and unexpected behaviour such as enhanced diffusion, breaking of self- averaging and 
emergence of Lifshitz tails |||| appear. 

All these latter phenomena were somehow believed to stem from the percolating lattice 
used to model the disorder. 

Here we shall investigate a different problem, where walks are partially and randomly 
absorbed at each site of a one-dimensional lattice thus leading to a non-conservation of the 
probability. From the physical point of view this model can mimic the partial absorption 
of a set of excitons wandering in mixed crystals. We shall show that despite the low- 
dimensionality of the model, a rich variety of features qualitatively similar to the ones 
observed in the trapping in percolating disorder can be found. The low dimensionality 
of the model however leaves open the possibility of a full analytical solution which would 
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extend the results presented here. 

Diffusion on a one-dimensional hierarchical lattice |TU] and multifractal characterization 
of the escape probability [TT] were previously considered in this context. 

Our aim here is then twofold. On one hand we shall show that most of the patholog- 
ical features of the model considered in Ref. ||[| are also present in this simplified low- 
dimensional version. On the other hand the results presented here can be regarded as a 
complement of earlier investigations f^,[n| where the total probability is conserved, thus 
providing a direct test on the effect of the non-conservation of the probability. 

The outline of the paper is as follows. In Sec. |I| we introduce the model and recall 



some well known general manipulations of disordered one- dimensional lattices. Sec. 1 1 1 
contains a numerical solution of the master equation. Sec. [IV] contains an effective medium 
approximation to the diffusion, which is shown to be inadequate. However in Sec. [V] an 
heuristic argument patterned after the Grassberger-Procaccia's similar reasoning for the 
Donsker-Varadhan case |7|], provides an intuitive explanation of the numerical results. Sec. 
|VT| contains a detailed numerical investigation of the influence of trapping on the spectrum 
and the localization properties. Finally in Sec. |V11] some conclusions are drawn up. 



II. THE MODEL 

Consider a particle moving on a one-dimensional lattice with random barriers and random 
trapping probability on each site. The master equation reads: 

(t) + ^x,x+l-Pxo,x+l (t) (1) 

for the probability P xo ,x{t) of being at site x at time t having started from site xq at the 
initial time t = 0. In this notation w X)V is the hopping probability from site y to site x, 
£x = w x -i iX + w x+ \ yX < 1 and 7 G [0, 1] is a parameter defining the sojourn probability which 
can be continuously tuned from 1 (no trapping) to (full trapping). Thus at each time step 
a particle at site x can move at sites x ± 1 with random probability w XjX ±i, stay at x with 
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sojourn probability 7(1 — ea,) and disappear with probability (l — j)(l — e x ). In the following 
we shall consider only the case of symmetric hopping probability (w x y = w VjX ) and infinite 
lattice. 

As it is well known when 7 = 1 this model can be mapped into a variety of other 
one-dimensional models. 

The (discrete) Laplace transform P XQtX (u) of the probability P Xo<x (t) satisfies an equation 
of motion which can be cast in the following Green equation for a Tight Binding (TB) 
hamiltonian: 

J2(E - Hi- )),,/;,, ,//•) = 5 XtX0 (2) 

y 

corresponding to the generalized TB hamiltonian 

74,3,(7) = Qx(l)5y, X - W X} y(5y iX -l + 5 y>x+1 ) (3) 

where we have defined ^(7) = 7e x + 1 — 7 and G XQ:X (E) = — -P C0 , a; (co , )| w =__B- Here the 
(discrete) Laplace transform is defined as: 

On the other hand eqn. (ffl) can be written in the following form: 

^ ,,(t + l) = £ T ^(7)iW*) (5) 

y 

where we have defined the transition matrix 



7(1-63;) ifx = y 

w x , y ]i\x-y\ = l (6) 
otherwise 



It is then easy to check that if 7 < 1 (and thus Y^ x T X}V < 1 for any x), then all the eigenvalues 
{A Q } of T are strictly less then unity, i.e. |A a | < 1 for any a, unlike the conserved case 
(7 = 1) where the fact that the maximum eigenvalue A m is non-degenerate and equal to 1 
irrespective of the (finite) size of the system, is ensured by the Perron-Frobenius' theorem 



p~3|1 . Equivalently this means that the eigenvalues {E a } of the Hamiltonian (|3|) are all strictly 
positive, i.e. E a > for any a. 

All the results presented in this paper were obtained using the discrete time equation 
([[]). The continuum time counterpart of eqn (|1|) would be, in the Coutinuos Time Random 
Walk (CTRW) formulation: 

dtP XQ ,x{t) = w x>x -iP XOiX -i(t) + w X)X+ \P XQ)X+ i(t) — e x P X0:X (t) — (3 x Px ,x(t) (7) 

where we have indicated with j3 x the absorption rate at site x and used the tildas to indicate 
the quantities which are rates. By balancing the gain, loss and sojourn terms to unity, it is 
easy to see that (3 X = (1 — 7)(1 — e x ). However it should be stressed that eq. (|7|) is not the 
continuum limit of eq. (|]). 

III. EXACT ENUMERATION 

We have considered a distribution of of the disorder given by: 

p(w) = 2 1 - a {l - a)w~ a 6{w)6(l/2 - w) (8) 

where a G (— oo, 1). By varying the parameter a we can pass from weak disorder (a < 0) 
to strong disorder (a — ► 1). The case a = corresponds to a uniform distribution and it is 
marginal in the sense that the inverse first moment is logarithmically divergent. We shall 
indicate with a overbar the average over the disorder (§) which is assumed to be quenched. 
The simplest quantity that one would like to compute is the mean-square displacements: 



>,,, . ^0 ) 2 Pxq, X (t) 



W ExPx ,x(t) ^ 

where it should be noted that the denominator must be included since the total probability 
is not conserved. The disorder average of the latter is called survival probability P s (t) and it 
is a decreasing function of time having chosen the initial conditions P s (0) = 1. The return 
probability 
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is another interesting quantity to look and it is again different from the survival return 
probability Pg(t) = P Xo ,x {t) in view of the non-conservation of the probability. A suitable 
normalization procedure || is to be used in order to avoid a quick underflow of the survival 
probability. 

We numerically solved the master equation exactly up to t = 1000 for three representative 
values of the strength of the disorder, namely a = —0.5 (Weak disorder), a = (Marginal 
disorder) and a = 0.5 (Strong disorder). The lattice size was chosen sufficiently large (up 
to iV = 2 18 ) to avoid finite size effect due to the boundaries. Our estimates are based on an 
average of 3 samples of 1000 different configurations each. Errors are statistical. 

The case 7 = 1 was previously studied in Ref. || and our numerical results are in perfect 
agreement with the one reported there. Indeed we find: 

< x \t) > = t 2v (11) 

where v is the correlation length exponent, and 

P ( t ) = t-W (12) 

where d s is the return probability exponent. Our estimates are 1v = 0.99 ± 0.01 and 
0.68 ± 0.01 for a = 0, 0.5 respectively which compares well with the expected values 1 and 
2/3. 

The case 7 = 0, corresponding to a full trapping, appears to be completely trapping- 
dominated. Indeed here the sample-to-sample fluctuations are enormous and they completely 
rule the diffusion. In order to have a feeling for this, we plotted in Fig. p] the mean- 
square displacement for a single configuration obtained with the same initializing seed, in 
the case 7 = 1 and 7 = for various strengths of the disorder. The staircase behaviour 
which can be observed in the case 7 = is a consequence of the non-conservation of the 
probability. A similar effect occurs in the return probability. This rather peculiar feature 



was already observed in the model of Ref . || . Here however it is noteworthy the extremely 
weak dependence of the behaviour case from the strength of the disorder. 

Upon disorder average we find, in the 7 = case, that the first moment (< x(t) >) is 
zero as expected and the second moment < x 2 (t) > and the return probability Po(t) follow 
the behaviour flTTD and flT2|) respectively. 

Our best estimate for the exponents v and d s are reported in Table |I[ As mentioned 
above, the very similar result for the three case of a is somewhat surprising, but it is in 
agreement with the aforementioned interpretation. Although the three values are not all 
within the numerical errors, we cannot rule out the possibility that the universality class 
be independent on the disorder choice. We also checked that all the intermediate cases 
< 7 < 1 behave as in the 7 = case after a transient time which depends on the value of 
7. It should be stressed that d s bears the meaning of spectral dimension || only in the case 
7 = 1. 

Although we monitored the behaviour of the survival probability and found that it decays 
as a stretched exponential (i.e. with an argument for the exponential which is a power of the 
time less than one), we also found that the numerical value of the exponent is not very easy 
to pin down. This was to be expected: a similar feature occurs in the Donsker-Varadhan 



problem [|14 



In view of these results, one expects both the survival probability and the survival re- 
turn probability to be non- self- averaging quantities. We investigated the full probability 
distributions of both quantities. We find (see Fig. |2|) that they rather accurately follow a 
log-normal distribution: 

V[X(t)) = —^eM ~ {hl *- Xt)2 ] (13) 

where X(t) = P X0 ,x {t)iYl,xPx Q .x{t)- Here Xt and o 2 are the mean and the variance of the 
distribution respectively. If asymptotically (t » 1) it happens that of >> X t then self- 
averaging is broken (see discussion in Ref. ||). Both these quantities can be computed 
directly as first (X t ) and second (a 2 ) moments of the distribution but also indirectly by 



fitting the evolution of the log-normal function at various values of t. We find the following 
behaviour for t » 1: 

o\ ~ t 2x (14) 

for the survival probability (and this defines the "free-energy" exponent x) an d a similar 
behaviour for the survival return probability (which defines the analog exponent Xo)- Our 
best estimates for these exponents is reported in Table |I[ On the other hand we find that 
asymptotically \ t ~ t both for the survival and survival return probability. Since in all cases 
2x > 1 the self-averaging property is broken ||. 

As a final remark we computed higher moments (< x A (t) >, < x 6 (t) >, etc) in the 
attempt to find signatures of multifractality as suggested by the results of Ref. [|ll[]. We 
found that all the moments were related to the second one irrespectively of the strength of 
the disorder, thus ruling out the possibility of multifractal behaviour in the size of the walk. 



IV. EFFECTIVE MEDIUM APPROXIMATION 

In this section we will tackle the problem of solving eq ([!]) by using an Effective Medium 
Approximation (EMA). Although this approximation is known to fail in many situations ||, 
it nevertheless gives extremely accurate results for the problem described by ([I]) with 7 = 1 
||. Here we shall carry out the successful recipe given in |J to the other extreme, namely 
the case of full trapping (7 = 0). 

The equation of motion for the Laplace transform P XOjX (u) can be written in the form 
(taking xo = for concreteness) 

Go, x (u) = /i x -i(uj)Qo tX -i(uj) + /i x (u)Q x+1 (u) + 5 X}0 (15) 

where Go tX {u)) = (l+u)P 0tX (u>) and where we have defined fi x (fi x -i) = w XjX+ i(w x -i iX ) / 
For x > and x < one can introduce the following new fields 

= ~~n 77a — ^ W = ~r 7777~ ^ 



respectively. In this way the problem is reduced to first-order and can be solved by continued 
fractions. For x > one finds the following recursions: 

A similar equation can be found for x < 0. It should be noted that the support for the 
have to be constrained in such a way that eq (|T7D be sensible. From eq. flip]) for the case 
x = and using eq. (|17D and the corresponding for x < one easily finds that 

Go *( u ) = 1 ^+77^ AT < \ ^ 
1 - ^{u) - 



The disorder average of eq. ([Tj|) thus gives: 
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0o,oM = Jd<P + J #_II+(0 + ) nj(0_) 9(1 - 0+ - 0_) (19) 

where the distributions 11+ (0+) and EE~(0_) are given by 

2 

n+(0 + ) = J dfi P {n) J d<p> + n+(0' + ) 5(<p + - (20) 

and similarly for H£ (</>-). Now comes the EMA approximation: 

n+(0 + ) = n-(0_) = 5(0-0 e H) (21) 

where <p e (uj) is an effective field to be found self-consistent ly. It is immediate to show that 
the self-consistency equation decouples and its solution is given by: 



M") = — v 2 (22) 

{^ a {ijj) = A*o/(l + &)) provided that the second moment of the distribution 

2 = I d}^ V 2 Pit*) / 9 oN 

Ma ~ /dA*p(A*) 1 } 

exists. As a consequence: 



SojoW = i (24) 



Due to conditions imposed on the transition rates, we have that \x a < 1/2. If equality 
holds then the result is identical to the one obtained in the absence of disorder. If however 
Ha < 1/2 then by antitransforming back to direct domain (in the continuum approximation), 
it is easy to show that the survival return probability behaves as: 

P°(t) = e~%{2n a t) (25) 



where Iq{z) is the 0-th order Bessel function. Upon asymptotic expansion |L5[ the leading 
behaviour for t » 1 is 

-(l-2Ha)t 1 1 

The result is qualitatively different from the numerical indication. This was to be expected in 
such [H EMA may not be capable of capturing the nuances of systems where large sample- 
to-sample fluctuations are present. 



V. AN HEURISTIC ARGUMENT 

In this section we shall present an intuitive argument yielding a stretched exponential 
for the survival probability in qualitative agreement with our numerical findings |16|| . The 
argument is based on a similar one given by Grassberger and Procaccia [f|J for the Donsker- 
Varadhan problem, and provides a plausible explanation of the irrelevance of the choice of 
the distribution for the hopping probability p(w). 

The main idea is that the maximum contribution to the survival probability is coming 
from rare regions where the hopping probabilities are all very close to 1/2, say in the interval 
[1/2 — e, 1/2] (with < e << 1/2) for the sake of the argument. If s is the number of sites 
of one of these regions, the typical number of steps r necessary to explore the region is 
r ~ s 1 ^ . If there were no leaking of probability the decay of the survival probability in 
this region P(t, s) would be exponential. However, during the time r, there is a loss of 
probability of order e T . The probability V(s) of finding such a region is (2e) s . Therefore the 
survival probability is expected to have the form: 
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P s (t) = / dsV(s)P(t,s) ~ / ds exp[~Ms)] (27) 

JO JO 

where for e << 1 we have that 

f t (s) = -L + \s^ + \s (28) 

where A = | lne|. The integral can be carried out using the steepest descent method which 
yields the following equation 

t = \ s 2 /u + \vs^ (29) 
for the saddle point sq. Since u < 1, at the leading order for t » 1, the solution is: 

s ~ (l/\) u/2 t u/2 (30) 
Therefore, upon substitution in eq. (|27l) we find at the leading order in t » 1: 

P s (t) ~ exp(-C v / t) (31) 

where C is a constant, which is different from the Donsker-Varadhan case where the argument 
of the exponential is ~ — t 1 ^ 3 . 

This argument seems to suggest the independence of the strength of the disorder which is 
irrelevant in the regions in the neighbourhood of 1/2. Clearly the argument can be extended 
to higher dimensions and self-similar lattices. 



VI. SPECTRAL ANALYSIS 

A. Density of States 

We shall here study the spectral properties of the TB equation associated to (fj): 

^( 7 )^(£) - Wavc-lVfe-l^) - Wx,x+lA+l(E) = E4j x (E) (32) 

by studying the density of states and the localization properties. For 7 = 1 the result is well 
known [T7]]. The Hamiltonian is positive definite as one can immediately check by defining 
the following creation 
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(o)^p) x = y/w XtX -ii/j x -i - ^/w XtX+1 ip x (33) 

and destruction 

(aip) x = yjw x+hx ip x+1 - y/w x , x+1 ij) x (34) 



operators and noting that 7^(1) = o)a [18|. Then for any 7 < 1 we have that, using < ip\(p > 
to indicate the usual scalar product: 

<^|H(7)|^> > < iJ)\H(i)\ip > >0 (35) 

and thus the Hamiltonian 7^(7) is strictly positive for any 7 < 1. 

A well known and efficient way to compute the spectrum of a one-dimensional disordered 
model is by means of the node- counting theorem |I7j]. By defining 



U x {E)=w x . ltX (36) 
Wx-i\E) 



the following recursion can be easily found from (|32|): 

up' 

U X+1 {E) = mi) ~ E) - (37) 

By calculation the number of times the quantity U X (E) changes sign, the integrated density 
of states 

M(E) = / dE' N(E') (38) 

(rather that the density of states N(E)) can be easily computed. The validity of the nu- 
merical procedure can be tested against the periodic case where the exact result 

1 1 

M(E) = arcsin(l - E) (39) 

2 7T 

with E G [0,2], can be easily derived. The case 7 = 1 was analytically solved by Stephen 
and Kariotis |i~2| . They found in the E — > limit: 



i?2-a for < a < 1 



M(E) = { (40) 

E 1 l 2 /J\\n.\fE\ fora = 
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Our numerical results for a = 0.5 and a = reproduce very well the theoretical prediction. 
Indeed we find 0.34 ± 0.01 and 0.45 ± 0.01 for the two cases where the theoretical values are 
1/3 and 1/2 respectively. The a = case is less precise due to the logarithmic corrections. 
Our results are based on an average of 5 different configurations with iV = 2 18 number of 
sites. We found that the the results were changing very little from one configuration to 
another (this is tantamount to say that N(E) is self-averaging) and thus this average was 
sufficient for our purposes. 

The behaviour for 7 = is qualitatively different (see FigJJ). We find for E — > (see 
Fig. |j) a Lifshits tail singularity: 

M(E) = exp(-l/E) (41) 



for all three cases a = —0.5, 0, 0.5. As expected on the basis of the result of Sec. |T|, all 
three values of the disorder give the same quantitative results. Fluctuations in the tail region 
grow as the disorder becomes stronger. We have been unable to give a theoretical derivation 
for this behaviour. It should be noted that the EMA of Sec. [TV] yields: 

N(E) = - - 1 = (42) 
V ; 7T V2£ - E 2 - c V ; 

for E G [1 — 2ji ai 1 + 2/i a ] and otherwise with c = 1 — 4/x^. With c = this is just 

the well known result for a periodic one- dimensional lattice. Again the EMA fails to give 

qualitatively correct predictions. 



B. Lyapunov exponent 



Eq. (|32|) can be cast in a transfer matrix form 



m x+1 {E)=T x {E)m x {E) 



(43) 



where we have defined the matrices: 



T X (E) 







1 \ 



Wx,x-l flx{l)-E 

\ Wx,x+1 / 



(44) 
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and 



ME) 



(45) 



/ 



Upon iteration eq. (^) becomes, if N is the number of sites: 

N-l 

9 N (E) = ( J] T X (E))V (E) 



(46) 



x=0 



where we have taken the initial vector as: 

ME) = 

The Lyapunov exponent is then defined as: 



V 1 / 



(47) 



7 (£) = lim ±hx\\* N (E)\\ 



(48) 



Again we tested our algorithm in the case 7 = 1 where an exact result [12] gives for E — > 0: 



7(E) 



E 2 -<* for < a < 1 



(49) 



E^ 2 /^\hiy/E\ fora = 

Our results are 0.34±0.01 (to be compared with 1/3) and 0.53±0.02 (to be compared with 
1/2). For E — the Lyapunov exponent is zero as it should since it corresponds to the 
stationary state of the master equation (|J) |1| . 



In the 7 = case instead we found that 7(0) = c(a) that is the ground state is localized 
(see Fig. |^). The spectrum is symmetric around E = 1 and this stems from the fact that 
the matrix 7i(0) — 1 is traceless. Again all the intermediate cases (0 < 7 < 1) follow this 
behaviour with a shift of the center of the spectrum which depends on the value of a. 



C. Exact diagonalization 



As a cross checking for the results of Sec. |Vl A| and | VI B| we direct diagonalized the 
matrix in @ both for the case 7 = 1 and 7 = by using standard IMSL package routines 
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for small matrices (N = 100). Although a direct quantitative comparison with the previous 
results is out of discussion since much higher sizes would be necessary, these results provide 
a qualitative understanding of the main difference in the spectrum of (|3]) in the cases 7 = 1 
and 7 = 0. In Fig. |6| the spatial dependence of the eigenvector ip x {E) is displayed for 
7 = 0, 1 for an energy in the lower (a), central (b) and upper (c) part of the spectrum. The 
chosen value for the strength of the disorder was a = 0, but no qualitative change in the 
behaviour is found for other values of the disorder in the relevant case 7 = 0. Correctly the 
ground state E = for 7 = 1 is extended as we argued before. 

The 7 = ground state is very localized consistently with our previous findings from the 
Lyapunov exponent. 

Another general probe to test the localization properties is to compute the Inverse Parte- 
cipation Ratio (IPR) 

v(E) _ HE x \ME)\r m 
p{E) ~n E.I^W (50) 

where E and ip x (E) are the eigenvalues and eigenvectors of eq. and N is the number 
of sites on the lattice. In this notation, then a state E is extended when p(E) ~ 1 while is 
localized whenever p(E) ~ 1/N << 1. The results are shown in Fig. ^ and are consistent 
with the previous picture. 



VII. CONCLUSIONS 

In this paper we have presented a detailed investigation of the properties of a one- 
dimensional disorder model for diffusion where at each site a fraction of the initial particles 
could disappear with random probability depending on the random adjacent barriers. Our 
work combined numerical and analytical techniques. The main results of this investigations 
can be summarized as follows. 

1) We showed that the trapped case (7 < 1) is qualitatively different from the conserved 
(7 = 1) counterpart and that a simple mean-field type of approach is not able to 
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capture the large fluctuation introduced by the non-conservation of the probability. 

2) We found that the trapped case appears to be completely disorder-dominated, thus 
suggesting a super- universal behaviour independent on the disorder (which is not the 
case for the conserved problem). We provided an intuitive explanation based on an 
heuristic derivation of the survival probability of why this may occur. 

We believe that more analytical work would be necessary for a complete understanding 
of the particular features appearing in the model considered in this paper. Work in 
this direction is ongoing and will be reported in the future. 
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FIGURES 

FIG. 1. Log-Log plot of the mean square displacement < x 2 (t) > for a single configuration 
in the cases 7 = 1 (no traps) and 7 = (full trap) for representative values of the disorder 
(a = —0.5,0,0.5). The 7 = case appears to be independent of the strength of the disorder 
parameter a. 

FIG. 2. Distribution density probability for the survival probability for time steps 
t = 200, 400, 600, 800, 1000. The full lines correspond to the fitted log-normal distribution of the 
form (|l3|). Shown are the cases a = —0.5,0,0.5 corresponding to (a),(b) and (c) respectively. 

FIG. 3. Integrated density of states on the whole spectrum E £ [0, 2] for 7 = 1 (a) and 7 = 
(b) for the same strengths of disorder as above. 

FIG. 4. Plot of the logarithm of the integrated density of states lnM(.E') versus 1/E for E — > 
in the case 7 = (full trapping) for a = —0.5, 0, 0.5. The behaviour is qualitatively different from 
the case 7 = 1 where M(E) is a power law in E. 

FIG. 5. Same thing as in Fig. ||, for the Lyapunov exponent j(E). A non-zero value of j(E) 
indicates that the eigenvalue E is localized. 

FIG. 6. Spatial dependence of the eigenvectors *p x (E) corresponding to the lower (Em) (a), 
central (Eq) (b) and upper (E m ) (c) part of the spectrum for the case a = and 7 = 0, 1. 

FIG. 7. Inverse Partecipation Ratio (IPR) defined in (|50|) as function of the eigenvalue E for 
a = and 7 = 0, 1. Values close to 1 indicate that the state is extended, while values close to 
(~ 1 /N) mean that the corresponding eigenvalue E is localized. 
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TABLES 

TABLE I. Summary of the exponents v and d s for the 7 = case in three representative values 
of the strength of the disorder. 



a 2v d s /2 

-0.5 1.24 ±0.01 0.57 ±0.01 

1.25 ±0.01 0.59 ±0.01 

0.5 1.26 ±0.01 0.60 ±0.01 



TABLE II. Summary of the exponents x an d Xo defined in the text. The label (D) and (I) mean 
direct evaluation and from the log-normal distribution, respectively. The values corresponding to 
a dash were considered unreliable due to the presence of a strong curvature in the preasymptotic 
regimes. 



« X(D) Xo(D) x(I) Xoffl 

-0.5 0.70 ±0.01 0.72 ±0.01 0.72 ± 0.01 0.70 ± 0.01 

0.73 ±0.01 0.72 ±0.01 0.70 ± 0.01 

0.5 - 0.72 ±0.01 0.70 ±0.01 
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